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Method of processing seismic data acquired by means of 
multicoioponent sensors 
GENERAL FIELD 

5 The invention relates to techniques for processing 
. seismic datal acquired by means of multi component 
sensors. 

This invention is in particular applicable to 
10 acquisition by means of cables disposed on the bottom 
of the sea (so-called ^^OBC or ^^Ocean Bottom Cable" 
techniques) . 

Multi component geophones capable of working in any 
15 position whatsoever, in particular at the bottom of the 
sea, have recently been proposed. This '"omnitilt" probe 
technology has allowed new simplified cables 
(mechanical joints are no longer necessary) and allows 
acquisitions with a better seismic bandwidth. 

20 

However, the acquisition step does not make it possible 
to provide the true orientation of the geophones of the 
cable, although this information is indispensable for 
making it possible to process the data. 

25 

The invention proposes a processing which is intended 
to be implemented on raw data and which allows 
reorientation and calibration (intended to convert the 
measurements of various geophones into a common phase 
30 and amplitude response) . 

STATE OF THE ART 

Techniques consisting in isolating from the signal the 
35 data which correspond to the first arrival at the 
sensor and in determining on the basis of these data a 
filter intended to be applied to the raw data so as to 
correct them and to thus obtain the components of the 
signal on the expected axes have already been proposed. 
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A proposal to this effect has been described in the 
article : 

""Horizontal vector infidelity correction by general 
5 linear transformation" - Joe Dellinger et al. - SEG - 
9-14 September 2001. 

However, this • technique is not necessarily optimal 
since the coupling mechanism v/hich intervenes at the 
10 geophone level is not the same for the waves which 
correspond to a first arrival at the sensor and for the 
waves reflected or converted by the seismic horizons. 

PRESENTATION OF THE INVENTION 

15 

The invention proposes another approach which employs 
the true data window for numerically reconstructing 
geophones oriented along the desired axes. 

20 Implicitly, this approach compensates for the er:rors 
which are not related to the geophones themselves, but 
which are due to the fact that the coupling between the 
geophone and the waves to be recorded is different 
depending on whether it is necessary to make a vertical 

25 vibration movement rather than a horizontal movement 
(on account of gravity) . 

In the case of a cable, the coupling is furthermore 
different depending on whether the vibration movement 
30 is in the direction of the cable or transverse. 

Moreover, since deeper windows are subject to a lower 
S/N (signal-to-noise) ratio, processing which 
implements trace stacks is moreover used. 
35 ... 

The invention thus proposes, according to a first 
aspect, a method of processing seismic data acquired by 
means of a sensor having at least three geophone 
components, characterized in that estimators are 



determined which are combinations of these components 
making it possible to. isolate the various data 
depending . on whether they correspond to propagation 
with reflection or with conversion and in that, to 
5 determine a sensor reconstruction,r the operators to be 
applied . to the various components of the sensor are 
determined in such a way as to . minimize the deviation 
between reference data and data obtained by applying 
the estimators to the sensor reconstruction, the 
10 operators thus determined being applied to the data 
acquired. 

It is specified here that, in the remainder of the 
present text, the term geophone is understood to mean 
15 any velocity sensor and the term hydrophone any 
pressure sensor. 

Preferred, but non limiting aspects of the method 
according to the first aspect of the invention are the 

20 following: 

the sensor furthermore including a hydrophone, the 
reference data for reconstructing a vertical geophone 
are derived from the data acquired by the hydrophone; 

the reference data for reconstructing a vertical 

25 geophone without hydrophone . or for reconstructing 
horizontal geophones are derived from the application 
of the estimators to one of the geophones of the 
sensor; 

the orientation in the . horizontal plane of a 
30 geophone component is obtained by minimizing the 
estimator of the transverse reflection; 

the estimators are determined as a function of a 
model of isotropic propagation or including the 
azimuthal anisotropy. 

35 

According to another more general aspect, the invention 
proposes a method of processing seismic data acquired 
by means of a sensor having at least three geophone 
components, characterized in that estimators are 



determined which are combinations of these components 
making it possible to . isolate the various data 
depending on whether they correspond to propagation 
with reflection or with conversion. The estimators thus 
determined may find applications other than that 
forming, the subject of the method according to the 
first aspect of the invention. 

DESCRIPTION OF THE FIGURES 

Figure 1 is a diagrammatic representation giving 
the angular conventions used; 

- Figures 2 and 3 are flow charts giving the various 
steps of the processing respectively in one and the 
other of the two exemplary implementations described. 

DESCRIPTION OF ONE OR MORE MODES OF IMPLEMENTATION OF 
THE INVENTION 

First exemplary implementairioii : case> of an iso-bropic 
propagation model 

Under the assumption of a locally ID (one dimensional) 
geology in proximity, to the receivers, and assuming 
isotropic propagation of the earth, a given geophone, 
with an orientation (j) \|/, measures: 

nik = Rpp cos {\|/)6ppjc + (Rps cos (Gk - ^) + Rtrsv sin(0jc - 
<|)) )sin(\j/)5psk 

With : 

k : index for the shotpoint (from 1 to N) 
0k : azimuth of the shotpoint with respect 

to the abscissa axis 
-Rpp. : reflectivity PP 

5pp : dynamic correction PP {^'normal 

moveout" or NMO) 
Rps : isotropic radial reflectivity PS 
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15 



20 



Rtrsv : isotropic 

reflectivity PS 
Sps : dynamic correction 
move out'' or NMO) 



PS 



transverse 



(^'normal 



This- model allows the evaluation of the reflectivity 
parameters from the set of traces tr^ through simple 
processing of least squares comparison (ignoring \|/ to 
begin with) in the Fourier domain, leading to the 
10 following equations: 



^ N ^ 



wc{(p) 







' Svx ^ 


lips 




Shcx{q>) 


^Rtrsvj 







Sc2{(p) Scs(<p) 
Scs((p) Ss2{^) 

Scalar quantities : 



Sc(<|)) = Zk cos (Gic-(|)) Ss(<|>) = sin (ek-(|)) 
Sc2((|)) - Sic cos^CGk-cI)) Ss2((j)) = sin2{ek-<|)) 

Scs ((|)) = Ejc cos{0jc-<|)) sin{6k-(|>) 
N = Sc2 + Ss2 



Wavelet quantities : 

wc(<|)) = Sk cos (0k-(|))5psk Sppk"^ 

25 ws((j)) = Sk sin(ek-<|>)5psk 5ppk"^ 

Trace stack quantities for qeophone x : 

Svx = Zk Xk 5ppk"^ 
30 Shcx(<j)) = 2k cos (0k-<}>) Xk 5pSk"'^ 
Shsx((|)) == Sk sin (Gk-<1>) Xk Spsk"'*' 

The solution of this linear system gives: 
Aiso Rpp cos (\|/) = W Svx + (Scs ws-Ss2 wc)Shcx + (Scs 
35 wc-Sc2 ws)Shsx 

Aiso Rps sin(v|/) = Kis cos (<|)) - Kic sin(<|)) 
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Also Rtrsv sin(\|/) = Kis cos((j)) - Kic sin (4)) 
With: 

W=Sc2 Ss2-Scs2 
5 Aiso = WN + wciScsyvs - Ss2 wc) + wsiScswc -^clws) 

Kic = {Scsw - Ss2 w)s vx + (n Ss2 - ws ws )sbcx + (- N Scs + ws w )shsx 
Kis = (Scs wc - Sc2 ws)Svx + (- NScs + \vc w )Sfacx +.(n Sc2 - wc wjshsx 

This modelling allows evaluations taking account of the 
following properties: 
10 a. Rpp does not depend on (|), 

b. 1 Rps I ^+ 1 Rtrsv I ^ does not of course depend on (|> 
either, 

c. Aiso is in practice rapidly steady over time and 
can. be ignored for the calibration/orientation 

15 procedure, since it is common to all the geophones of 
one and the same receiver. 

Evaluations of dense shots 

20 Most of the OBC acquisitions are gleaned using a dense 
and regular grid of sources, which allows considerable 
simplification: 

Sc = Ss = 0, wc = ws=0 (symmetry of the sources with 
respect to the receivers) 
25 Scs = 0 Sc2 = Ss2 = N/2 (isotropic source distribution) 

Next, the exact solution can be obtained through the 
approximation: 
N Rpp cos (\)/) = Svx 
30 N Rps sin(\|/) = 2 Shcx(<|)) 
N Rtrsv sin(\|/) = 2 Shsx((|)) 

This approximation leads to very simple calculations, 
not involving any wavelets, and can be applied 
35 immediately. 



Orientation of the geophones 
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Since Rtrsv does not exist physically, the minimization 
of the energy of Rtrsv leads to a trigonometric 
equation which gives the true orientation 4>geo( + ^ ^) • 

5 tan (2 (|)geo)= 2 iZ^o, Kis,)/{Y.Yii^ ^j;^Kis,') 

( (Emax-Emin) / (Emax+Emin) ) gives ai check on the 

quality of the reorientation. 

10 Moreover, if one wishes to find the orientation 
according to the first arrivals, it is possible to 
correct the said first arrivals so as to set them to 
one and the same arrival time, then to simplify kic and 
kis by replacing the wavelets wc and ws by the scalars 

15 so and ss, by considering that the waves recorded 
horizontally are in fact the projection of the radial 
wave P, present on all the geophones since it is 
oblique. . 

20 Geophone vertical composite calibration : 

With the' geophones gi, gz, ga^ we construct a vertical 
composite geophone v, v = opi*gi + op2*g2 + op3*g3 (or 
comprising additional similar terms in the case where 
25 extra geophones are present in the receiver) where opi, 
op2/ op3 are the filters of finite length and opu*gu 
represents the convolution of geophone gu with filter 
opu. 

30 such that: 

El = |XH - XVI^ - |Kic(v) 1^ + IKis(v) 1^ 

The energy of the difference between XH (hydrophone 
35 after application of the geophone phantom, or 
cross-ghost hydrophone) and XV .(the vertical composite 
geophone after application of the hydrophone phantom or 
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cross-ghost geophone) , (see for example in this regard 
the Applicant's Patent Application FR 2 743 896). 

E2 = |Rps(v)|^+ IRtrsv(v) |- 
5 horizontal energy of the vertical composite, 

E = XEI + (1-X)E2 is a quadratic form in the 
coefficients of the filters and can be reduced to the 
minimum, thus giving a linear system to be solved. {X 
10 is a matching parameter, 0<=X<=1, which favours either 
a greater adjustment to the reference hydrophone or a 
greater minimization of the shear energy) . 

In the case of terrestrial data, that is to say if 
15 there is no hydrophone available, it is possible to 
choose one of the geophones as reference and to replace 
the hydrophone by Rpp(gref)- 

Calibration in a horizontal arbitrary direction 

20 

With g = opi*gi + op2*g2 + op3*g3 and <pg an arbitrary 

direction. 

We define: 

El = |Rps(g, 9g) - Rps(ref, cpref)|^, 
25 as being the energy of the difference between the 
evaluation of Rps of the arbitrary composite geophone 
and the evaluation of Rps of a reference geophone (in 
general the geophone oriented in the direction of the 
cable) . 

30 

E2 = |Rpp(g)p + |Rtrsv(g, 9g)|^ (the nonradial energy) 

E = XEI + (1 - X)E2 allows the derivation of a 
composite horizontal geophone in the desired direction, 
35 having the- .same frequency response as the reference 
geophone, and with a minimum PP contamination. 
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Considering the cases (pref = 0 9ref = 7r/2, it is 

possible by simple trigonometric combination to 
generate the radial and transverse projections. 

5 Second exemplary implemen-ba-tion : case of anisotropic 
azimuthal propagation modelling 

Let a be the direction of the natural fast propagation 
axis. The modelling of the measurement of the geophone 
10 becomes (using one or other of Rpsl and Rps2 the two 
images along the natural directions, or Rps and 6Rps 
defined by Rpsl = Rps + 5Rps, Rps2 = Rps - SRps) : 

mjc = Rpp cos(vj/)6ppk + (Rpsl cos (cp-a) cos (9):-a) + Rps2 
15 sin (cp-a) sin (9k-a) ) sin (v|/) 5psk 

nik = Rpp cos (\t/) Sppk + (Rps cos(9k-q>) + 5Rps cos (6^+9- 
2a) ) sin (\|/) 5ps*' 



20 



giving the normal equations 



M = 



wc 



N 



wc 



M. 



((p) wc(2a-<p) 

Sc2(q)) Sc2(a)-Nsin^(a-{p) 
ws(2a-cpj Sc2(a)-Nsin^(a-((») Sc2(2a-9) 

^ svx ^ 

Shcx(q)) 



Rpp 
Rps 

SRps 



V 



Shcx(2a-q)) 



The solution of this .linear system gives: 

Also Rpp cos (y) = unchanged 
25 Aaniso Rps sin(\|/) = (Kac cos (2a-(p) + Kas sin ( 2a- 
(p) ) sin (2 (a-cp) ) 

Aaniso Rtrsv sin(\j/) = (-Kac cos (cp) -Kas sin ((p) ) sin (2 (a- 
(?) ) 

30 With: 
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Aaniso = sin^ (2 (a-9) ) Also 

Kac = (Sc2 w-Scswc ) Svx + (N Scs-wc m'j ) Shcx - (N Sc2- 
wcwc ) Shsx 

5 Kas = (-Ss2 wc+Scswj)Svx + (N Ss2-ws vw" ) Shcx - (N Scs- 
ws wc ) Shsx 

Vertical calibration of compositie geophone: 

10 The isotropic process remains applicable with the 
change 

E2 = |Kac(g)p + |Kas (g)p 

Horizontal arbitrary calibration or rows /columns of 
15 composite geophones 

The observation of SRps over the data field makes it 
. possible to diagnose the presence {or otherwise) of 
significant azimuthal -anisotropy. (the quantity 
20 sin^ (2 (a-cp) ) SRps does not require a knowledge of a , f or 
its calculation) . 

The isotropic process remains applicable with the 
changes 

25 E2 = |5Rps{v)|2 and E = X (E1+E2 ) + ( 1-A.) E3 . 

When a is not generally known, a scan over a range of 
7i/2 is implemented, using the value of a which 
minimizes Emini/EO. 
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CLAIMS 

1'. Method of processing seismic data acquired by 
5 means of a sensor having at least three geophone 
components, characterized in that estimators are 
determined which are combinations of these components 
making it possible to isolate the various data 
depending on whether they correspond to propagation 

10 with reflection or with conversion and in that/ to 
determine a sensor reconstruction, the operators to be 
applied to the various components of the sensor are 
determined in such a way as to minimize the deviation 
between reference data and data obtained by applying 

15 the estimators to the sensor rieconstruction, the 
operators thus determined being applied to the data 
acquired - 

2. Method according to Claim 1, in which, the sensor 
20 furthermore including a hydrophone, the reference 'data 

for reconstructing a vertical geophone are derived from 
the data acquired by the hydrophone. 

3. Method according to Claim 1, . in which the 
25 reference data for reconstructing a vertical geophone 

without hydrophone or for reconstructing horizontal 
geophones are derived from the application of the 
estimators to one of the geophones of the sensor. 

30 4. Method according to Claim 1, characterized in that 
the orientation in the horizontal plane of a geophone 
component is obtained by minimizing the estimator of 
the transverse reflection. 

35 5. Method "according to one of the preceding claims, 
characterized in that the estimators are determined as 
a function of a model of isotropic propagation or 
including the azimuthal anisotropy. 
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6. Method of processing seismic data acquired by 
means of a sensor having at least three geophone 
components, characterized in that estimators are 
5 determined which are combinations of these components 
making it possible to isolate the various data 
depending on whether they correspond to propagation 
with reflection or with conversion. 
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ABSTRACT 



The invention relates to a method of processing seismic data acquired by means of a 
sensor having at least three geophone components, characterized in that estimators are 
determined which are combinations of these components making it possible to isolate the 
various data depending on whether they correspond to propagation with reflection or with 
conversion. The estimators find application in particular for determining a sensor 
reconstruction according to which the operators to be applied to the various components of 
the sensor are determined in such a way as to minimize the deviation between reference data 
and data obtained by appljdng the estimators to the sensor reconstruction, the operators thus 
determined being applied to the data acquired. 
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Oriented horizontal geophone (pg 
g = opi'gt +0p2*g2+0p3*g3 



Choice 



i 



of gref I 



Yes 



Orientation of gref known? 



Minimization of Rtrsv(gref)|* 
with respect to <Pref 



El = |Rps(g, <Pg) • Rps(gref, <Pref)h 
E2 = |Rpp(g)p -».|Rtrsv(g, (pg)|^ 



Minimization of E = XE1 + (1 -;i)E2! to detennine op^ , op 2, 0P3 



FIG.2 
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Vertical geophone g = op^ +op2*g2 +op3*g3 



Hydrophone known? 



No 




XH = Slack Hydrophone* phantom_geophone 
El = I R pp (g* phantom.hydrophone ) - XH| ^ 



Choice of gref 
G = Rpp{gref) 
El - |Rpp(g)-G|V 



E2 = jRps(g)|^ + lRtrsv (g)|^ 



Minimization of E = X.E1 + (1-X)E2 to determine op^, op2, 0P3 



FIG.3 



